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Abstract: 

In  this  paper  we  introduce  a  new  definition  of  super  convergence  -  the  17%- 
sup trconv erg trice ,  which  generalizes  the  classical  idea  of  superconvergence  to  gen¬ 
eral  meshes.  We  show  that  this  new  definition  can  be  employed  to  determine 
the  regions  of  least-error  in  any  element  in  the  interior  of  any  grid  by  using  a 
computer-based  approach.  We  present  numerical  results  for  the  standard  displace¬ 
ment  finite  element  method  for  the  scalar  equation  of  orthotropic  heat-conduction, 
for  meshes  of  conforming  triangles  of  degree  p,  1  <  p  <  5,  and  elements  in  the 
interior  of  the  mesh.  The  results  demonstrate  that,  unlike  classical  superconver¬ 
gence,  17%-superconvergence  is  applicable  to  the  complex  grids  which  are  employed 
in  practical  engineering  computations. 
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1  Introduction. 


Let  be  a  one  parameter  sequence  of  finite  element  solutions  of  a  problem 
which  are  computed  using  a  sequence  of  meshes  T  =  {Th}  and  let  u  denote  the 
exact  solution.  Let  us  assume  that  we  are  interested  in  the  values  of  the  solution  or 
its  derivatives  or  linear  combinations  of  these  quantities  i.e.  in  the  linear  functional 
F(v)(m).  Let  us  assume  that  for  every  element  r  of  the  mesh  Th  a  special  point  e, 
which  depends  on  the  geometry  (but  not  the  size)  of  the  element,  is  given.  Then 
denoting 


(1.1)  ¥(tt  -  tih)  :=  max  \F(u  -  uh)(*)| 

we  are  interested  in  the  values  of  relative  error  in  F(u)  at 


(1.2) 


0(*;F;tt,ttfc,fc,r)  := 


0 


if¥(«-t*fc)/0  ; 
if  ¥(ti  -  tx  J  =  0  . 


If  the  point  m  is  such  that 


(1.3) 


e(*;  t)  <  jq0  *  —  ° 


then  m  will  be  called  a  tt-q%-  superconvergence  point  relative  to  the  exact  solution 
u  and  the  family  of  meshes  T.  Consider  now  a  family  U  of  solutions;  the  point  m 
will  be  a  superconvergence  point  relative  to  the  family  of  meshes  T  if  it  is 

tt-q%-superconvergence  point  for  every  u  G  U.  Obviously  ©(»; F;u,uh,h,r)  <  1 
and  thus  all  points  in  every  element  r  are  100%-super convergence  points.  Note 
that  if  there  exists  a  point  e  in  the  element  r  for  which  if  =  0  in  (1.3),  i.e. 


(1.4)  0(*;F;u,tth,fc,r) — >0  as  h — >0 

for  a  particular  u  (reap,  for  every  u  G  U)  then  m  is  a  u-super convergence  (resp. 
W-superconvergence)  point  tn  the  classical  sense.  The  superconvergence  of  this 
type  was  studied  by  many  investigators;  see  for  example  [1-28]  and  the  citations 
in  these  papers.  It  should  be  noted  that  the  q%-8uperconvergence  has  not  been 
analyzed  in  the  literature.  In  this  paper  we  will  study  the  *7 %-super convergence, 
when  F(v)  is  a  linear  combination  of  derivatives  of  u,  for  the  model-problem  of 
orthotropic  heat-conduction  for  meshes  of  triangles  of  degree  p,  1  <  p  <  5.  In 


[29]  we  will  address  17  %-super  convergence  for  meshes  of  quadrilaterals  with  local 
refinements. 

In  [28]  we  studied  classical  superconvergence  by  a  computer-based  approach. 
We  observed  that  the  superconvergence  points  are  very  sensitive  to  the  geometry  of 
the  mesh,  the  solution-type  and  the  coefficients  of  the  differential  operator.  Hence 
the  superconvergence  points  may  not  exist  (in  the  classical  sense)  for  the  complex 
grids  employed  in  practical  computations.  In  contrast,  the  77  %-super  convergence 
points  will  always  exist  for  a  sufficiently  high  17  and  can  be  defined  for  entire  classes 
of  grid  geometries,  solution-types  and  coefficients  of  the  differential  operator.  Thus 
the  ij%- superconvergence  points  are  robust  and  can  be  used  to  sample  the  solution 
m  practical  engineering  computations. 

Following  this  Introduction  we  outline  the  model  problem  of  orthotropic  heat- 
conduction  and  we  define  various  quantities  associated  with  the  geometry  of  the 
error.  We  outline  the  theoretical  setting  of  this  study  and  give  a  computer-based 
approach  to  obtain  the  q%-superconvergence  points.  We  then  give  various  exam¬ 
ples  of  the  77  %-super  convergence  points  for  various  grid  geometries,  solution-types 
and  coefficients  of  orthotropy  and  demonstrate  how  to  obtain  robust  sampling 
points  for  the  practical  grids. 

2  Preliminaries. 

2.1  The  model  problem 

Let  us  consider  the  model  problem  of  heat-conduction  in  orthotropic  medium 
with  mixed  boundary  conditions.  Let  fi  C  RJ  be  a  bounded  polygonal  domain  and 
let  its  boundary  dtl  be  split  into  two  parts  I'd  and  IV  (I'd  has  positive  length). 
Let  u  be  the  solution  of  the  problem 


(2.1a) 

£<*)!=  =/> 

(2.16) 

u  =  0 ,  on  I'd  » 

(2.1c) 

a 

2  =  9 »  on  TN  . 

)«=i 


Here  Ku,  k,l  =  1,  2  are  the  entries  of  the  thermal-conductivity  matrix  which 
satisfy  the  conditions 


(2.2«)  Ku^Ka,,  k,i  =  1,2, 

(2.26)  0  <  K^tl  +  i l )  <  £  KuUti  <  K^Ul  +  g)  V(  =  (6,6)  €  Ra , 

k4=i 

where  JL..  K —  are  the  principal  thermal  conductivities; 

(2.3)  «t(«)  ■='t,Xup-  ,  *  =  1,2, 

£=1 

axe  the  components  of  the  flux  (heat-flu*)]  n*,  6  =  1,  2  are  the  components  of  the 
unit  outer  normal  of  0fl;  /  €  L*(0),  g  6  £3(IV)  are  given  data. 

Let  H$B(tt)  :=  |u  €  H^ft)  |  u  =  0  on  I'd}.  Then  the  above  problem  may  be 
put  in  the  variational  form:  Find  u  G  lff0(Q)  such  that 

(2.4)  4»(«h«)  =  La(v)  Vug  j£,(Q) 
where 


(2-5)  BnM:=  L£Kuwt£;'  L°w:=Jjv+LSv- 

Let  T  :=  {Ta}  be  a  regular-family  of  meshes  of  triangles  with  straight  edges. 
(It  is  assumed  that  for  any  triangles  r;,  r,-  G  Ta,  the  intersection  D  r,-  is  either 
empty,  a  vertex  or  a  common  edge,  and  that  the  minimal  angle  of  all  the  triangles 
is  bounded  below  by  a  positive  constant,  the  same  for  all  meshes.)  The  meshes  Ta 
are  not  assumed  to  be  quasiuniform.  We  introduce  the  conforming  finite-element 
spaces 

(2.6)  SJ:=  {.€  C*(n)|  «UE  ?*(«),  *  =  1 . «(2\)}, 

where  V9(r)  denotes  the  space  of  polynomials  of  degree  p  over  the  triangle  r; 
M(Th)  is  the  number  of  elements  in  7*.  The  finite  element  solution  uh  of  the 
model  problem  satisfies:  Find  uh  G  5£rc  :=  fl  H£D((l)  such  that 

(2-7)  Ba(uh,vh)  =  La(vh)  V  vh  G  ^,rD 

We  let  eh  :=  «  —  i*fc  denote  the  error  in  the  finite  element  solution. 
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2.2  Definition  of  t?%-superconvergence  quantities 


Let  Th  6  T  be  a  finite  element  grid,  r  €  Th  be  any  element  and  K  be  the 
thermal-conductivity  matrix  expressed  in  a  given  fixed  coordinate  system.  Let 

F(u)  be  the  solution-quantity  of  interest,  for  example  we  may  have  F(u)  =  — — 

du  .  du  du  ,  .  du  du 

or  aT, "  *<*) = + = **  + *»  &;- e,c- We  B0W 

define  several  geometrical  quantities  associated  with  the  error  in  the  finite  element 
solution  in  the  element  r. 


Given  q,  0  <  q  <  100  we  define  the  following: 


1.  rj%-contour  of  F(u)  in  the  element  r  €  Th  for  the  exact  solution  u  and  for 
material- orthotropy  K: 


(2-8)  <$„,(«;  if  ir,Tk)  :=  {«  €  r  ]  6(*;F;K;*,uh,h,r)  =  ^  }  , 

Here  0(«;  F]K]u,uk,h,r)  is  the  relative  error  as  defined  in  (1.2);  note  that  K 
was  added  to  the  list  of  arguments  to  indicate  the  dependence  of  the  relative  error 
on  the  orthotropy. 

2.  ri%-band  of  F(u)  in  the  element  r  €  Th  for  the  exact  solution  u  and  for  material- 
orthotropy  K: 


(2.9)  B^u)(u]  K-,r,Th)  :=  j®  €  r  |  B(m-}F;K]u,uhih,T)  <  } 

S.  Superconvergence  points  of  F( u)  in  the  element  r  G  Th  for  the  class  of  exact 
solutions  U  and  for  material-orthotropy  K: 


(2.10)  XXAU-, K; r,  Tk)  :=  0  r, Tk) 

u£U 

4-  rj%-superconvergence  regions  of  F(u)  in  the  element  r  G  Th  for  the  class  of 
exact  solutions  U  and  for  material-orthotropy  K: 


(2.11)  *#„,(«;  K-,  r,  T*)  :=  f)  »>*)(«;  K<  *.  Tk) 

«€  U 

5.  Common  tj%-superconvergence  region  of  .F(tt)  in  the  element  r  G  TK  for  the 
class  of  exact  solutions  U,  the  class  of  materials  M  and  the  class  of  meshes  T : 
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(2.13) 


*5So(KiA4iT,T)=  n  n  *&)(«: Jr;T,rj 

TkeT  KSM 

Here  M  is  a  class  of  thermal  conductivity  matrices  referred  to  a  given  fixed  co¬ 
ordinate  system  and  r  is  understood  to  be  an  element  of  fixed  geometry  (but  not 
fixed  siae)  embedded  in  all  the  meshes  in  the  family  T. 

Remark  2.1.  In  some  cases  the  sets  defined  above  may  be  empty. 

Remark  2.2.  The  function  u  is  a  solution  (resp.  U  is  a  class  of  solutions)  of  (2.1) 
for  a  given  set  of  data  (resp.  for  given  classes  of  data).  When  /  =  0  we  say 
that  exact  solution  u  is  “harmonic”  (resp.  U  is  the  class  of  “harmonic”  solutions). 
(When  u  satisfies  Laplace’s  equation,  —  At*  =  0  (i.e.  when  K  =  J),  the  solution 
is  harmonic  in  the  classical  sense.)  The  majority  of  steady-state  computations  in 
engineering  are  done  to  approximate  “harmonic”  solutions;  thus  it  is  important  to 
study  superconvergence  for  this  class  of  solutions. 

One  is  interested  to  know  a-priori  the  ^%-super convergence  regions  for  the 
smallest  possible  values  of  q%,  for  which  they  exist,  for  classes  of  solutions  of 
interest  (In  practical  computations  in  plane  elasticity  and  heat-conduction  the 
class  of  solutions  of  interest  is  the  subclass  of  “harmonic”  solutions  with  a  finite 
number  of  algebraic  point  singularities  of  the  type  r“).  In  general,  if  Th  is  any 
grid,  it  is  impossible  to  predict  the  locations  of  the  q%-superconvergence  regions. 
Here  we  will  make  additional  assumptions  about  the  approximation  which  will 
enable  us  to  determine  a-priori  the  q%- super  convergence  regions  for  the  meshes 
which  occur  in  engineering  computations.  We  will  consider  that  the  problem  was 
solved  mriwg  a  sequence  of  meshes  {Th}  (not  necessarily  uniform  meshes)  which 
are  constructed  by  a  mesh-generator  or  an  adaptive  process.  We  will  address  the 
superconvergence  as  h  — ►  0  i.e.  we  are  interested  in  the  asymptotic  locations  of  the 
ti%-superconvergence  regions  in  the  limit  as  the  mesh-parameter  h  tends  to  zero. 

In  this  paper  we  will  determine  the  asymptotic  17  %-  super  con  vergence  regions 
for  elements  in  the  interior  of  the  mesh  and  smooth  solutions.  (The  solutions  can 
have  isolated  point-singularities  at  the  boundaries  but  they  are  analytic  in  the 
interior  of  the  domain.)  The  theoretical  results  will  be  given  for  a  special  class  of 
locally  periodic  meshes.  However  through  numerical  examples  we  will  demonstrate 
that  the  conclusions  of  the  theoretical  study  hold,  for  all  practical  purposes,  for 
the  complex  grids  which  are  used  in  the  practical  engineering  computations  (for 
example  the  grid  shown  in  Fig.  1). 

2.3  The  class  of  locally  periodic  meshes 

We  now  present  the  definition  of  a  special  class  of  locally  periodic  meshes.  Let 
us  consider  a  locally  periodic  meshes  (or  grid)  defined  as  follows.  Let  0  <  H  <  H° , 
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*5)  e  n, 

(2.13)  S(»0,fT)  :={•  =  (*„*,)!  i  =  1,2} 

and  assume  fT°  is  sufficiently  small  such  that  5(«°,  if0)  C  (1.  Further,  let  7  be  a 
set  of  multi-indices  (i,j) ,  e  0  and 


(2.14)  c(«W>,  k)  :=  S(«W>,  k)  C  S(»°,  B) ,  («,,•)  6  7 

be  the  set  of  the  h- cells  (or  cells)  which  cover  exactly  5(a°,  H)  i.e. 

(2.15a)  (J  c(*M\h)  =  §(m°,H) 

(2.156)  c(^\k)f]^\k)  =  0  for  (i„Ji)  ^ 

We  will  refer  to  S(m°,  H)  as  the  subdomain  of  periodicity  of  the  mesh  centered  at 
»°.  Denoting  by 


(2.16)  c  :=  5(0,1)  |(*i,ia)|  |*i|  <1,  |*a|  <  l} 

the  unit-  (master-)  cell  e,  the  A- cell  is  an  h-scaled  and  translated  master-cell. 

Let  T  be  a  triangular  mesh  on  the  master-cell  (the  master-mesh)  and  be 
the  mesh  on  d(m^\h)  which  is  the  scaled  and  translated  image  of  T.  We  will 
consider  the  family  T  of  locally  periodic  meshes.  Let  Th€T  and  Th(m°,  H )  be  the 
restriction  of  Th  on  S(m,H)  and  the  restriction  of  Th(m°,H)  on  c(»(’*J),h). 
We  assume  that  ,  (»,  j)  E  7  i.e.  Th(m°,H)  is  made  by  the  periodic 

repetition  of  the  h- scaled  master  mesh. 

The  type  of  meshes  under  consideration  is  depicted  in  Fig.  2a,  where  the  subdo¬ 
mains  5(ss°,  H)  D  S(m0,Hx)  are  covered  by  a  periodic  array  of  cells.  An  example 
is  shown  in  Fig.  2b  (where  the  master-mesh  in  the  master-cell  is  shown  in  Fig.  2c). 
Here  we  are  interested  in  the  superconvergence  regions  in  the  elements  belonging 
to  the  patch  shown  in  the  Figs.  2b,  2c.  Outside  the  square  5(*°,  H)  the  mesh  is 
arbitrary;  it  could  have  curved  elements,  etc.  The  values  of  Hx  depend  on  h  i.e. 
Hx  =  Hx(h).  More  precisely  we  will  assume  that  there  exist  Cx  and  C7  indpendent 
of  h  such  that 

(2.17)  CXH*  <h<  C7HX ,  a  >  1 
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3  The  theoretical  setting 

We  will  here  for  simplicity  address  only  the  case  of  Poisson’s  equation  (K  =  I 
in  (2.1)).  Nevertheless,  the  results  could  be  easily  generalized  to  the  general  setting 
for  orthotropic  medium. 

Given  a  function  it  and  the  multi-index  a  :=  (ai,a2)  we  define 


(3.1a) 


D°u  := 


0*?1  dxV 


,  \ct\  :=  +  a3 


(3.16)  (2?*it)(»)  :=  2  ^  »  k>0,  integer 


and  let  us  denote 


(3*2a)  |  |u|  15^^ ) 1 1«| lri(s(««jBr. )) 

(3*26)  I u  I  sw.Ht ) :=  llttIU-(S(»«^. )) 

(3.2c)  ll«lls(^jr,.)  :=  E  l^“«l 

|«|=i 


Let  Q  be  a  polynomial  of  degree  p  +  1  on  the  master-cell  c  and  let  T  be  the 
master-mesh.  Then  denote 

(3.3)  p:=  Q-  Q?T 

where  0^  is  the  interpolant  of  degree  p  of  the  function  Q  defined  over  the  master- 
mesh  f  (for  which  h  =  1).  Then  we  have 


(3.4) 


p(l»*a)  —  /*(  »  l*a I  ^  1 

p(*l,1)  =  p(*1,-l),  I*! I  <  1 


9 


because  Sf(c)  includes  all  polynomials  of  degree  pone.  Further  let 


(3.5)  ^pbr («)  :=  |u  €  Bl(c)  |  v  satisfies  (3.4)  j 

and 

(3.6)  S{m(c)  :=  {«  €  JSTJ„(c)  |  «  €  V,{r)  V  f  €  f  } 

Further  let  zp  €  5f  PER(c)  sue-  1 

(3.7o)  Bt{y,v)  (/>,$)  V  v  6  5T,pbr(c) 

and 

(3.76)  ^(p-i')=  0 

Note  that  the  function  V  exists  and  is  uniquely  determined  (we  will  compute  it 
numerically  in  the  examples).  Let  us  also  define  €  Bl{c)  by 

(3.8)  V>  :=  p  —  zp  =  Q  —  —  zp  . 

Below  we  outline  the  main  theorem  of  the  paper. 

We  will  make  the  following  assumptions  about  the  exact  solution  it: 

Assumption  I 

On 

(3.9)  |  D*u  |  <  K  <  oo,  0  <  |a|  <  p  +  2 

Remark  S.l.  Assumption  I  states  that  the  solution  is  locally  smooth  in  the  sub- 
domain  S(m°,  H)  i.e.  the  subdomain  should  be  sufficiently  far  from  boundaries, 
material-interfaces  and  points  where  the  data  are  rough. 

Assumption  II 

If  oa  :=  (D*ti)  (®°),  a  =  (aj,  a3) ,  )a|  <  p  + 1  then 


(3.10)  i?2=  £  a*  >  0 

W=p+i 

Farther,  we  outline  that  the  mesh  T( ft,  h)  is  such  that: 
jluitmption  III 
On  5(*°,fT1), 

(3-11)  IMI*>*>  <  C4' a. 

with  (}  >  (p  +  1)  —  e,  where  e  is  specified  in  the  theorem  below  and  where  C  is 
independent  of  T(ft,  A),  Hi,  bat  it  depends  on  K  and  R. 

Remark  S.S.  We  do  not  assume  that  u  is  smooth  in  ft  outside  of  5(as°,  H).  For 
example,  ft  can  have  a  boundary  with  reentrant  corners  (as  in  Fig.  2(a))  and 
hence  «  can  be  unsmooth  in  the  neighborhood  of  these  comers.  Nevertheless 
assumption  III  makes  an  implicit  requirement  on  the  (refinement  of  the)  mesh  in 
the  neighborhood  of  these  comers.  If  u  is  smooth  in  a  convex  ft  and  the  mesh  is 
quasi-uniform  then 

(3.12)  IM„  <C**-1Mr||2J'M«||„,  r  >  0 

and  hence  in  assumption  III  we  can  take  fi  <P+  1  arbitrary.  Assumption  m  is 
related  to  the  pollution  problem  which  is  analyzed  in  [38]. 

Remark  S.S.  Assumptions  II,  III  imply  that  the  principal  part  of  the  error  in 
S( ®°,  H7)  is  related  to  the  non-zero  (p  +  l)-derivatives  of  the  exact  solution  at  »°. 

Let  h))  be  the  function  rf),  defined  above,  scaled  and  trans¬ 

lated  onto  the  cell  c(a^,h)  of  the  mesh  in  S(m°,H)  i.e. 

(3.13)  4>dm):=  '!’(*),  §^(®) =  hP  ’  *  =  1’2’ 

where  m  €  c(*^,h) ,  a  =  i  (»  -  b^)  . 

h 

It  is  easy  to  see  that  i>h  can  be  periodically  extended  over  5(b°,JET1).  In  [28] 
we  have  proven  the  following  theorem  using  the  theory  of  interior  estimates  (see 

[31]-[37]): 

Theorem  1 .  Let  H7  <  Hx  <  H  <  H°  and  the  assumptions  I-III  and  (2.17)-(2.18) 
hold  with 


1 
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(3.14) 


6p  +  1 
6 p  ’ 


0  =  P+  1-6, 


6=1/ 


1 

6(6p  +  l) 


Then  for  any  a  6  S(m°,H7) 

(3.X5) 

with  |A|  <  1  and  C  independent  of  h. 

Remark  8.4.  Note  that  (and  thus  i/)h)  depends  on  the  polynomial  Q  which 
approximates  well  the  solution  in  5(as°,  H).  Theorem  1  states  that  for  small  h 
(and  hence  H7)  we  have  on  5( ss°,  H7) 


(3.X6)  ^(.)  *  g(.) 

or  more  precisely 

(3.X7) 

On  the  other  hand  it  is  easy  to  see  that 

(3.18)  I&U-J  >-Chr 

where  C  depends  on  the  constant  R  defined  in  (3.10).  Nevertheless  this  does  not 
mean  that 

(3.19)  l§rl.- c*'-  *  =  1,2 

The  case  that 


(3.20) 


1 

<r0  >  0  for  either  *  =  1  or  2  I 
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is  very  exceptional.  Hence  we  will  assume  that  (3.20)  does  not  happen  i.e.  we  will 
assume  that  there  exists  a  constant  £  >  0  independent  of  h  such  that 


(3.21) 


This  assumption  could  formally  be  achieved  by  imposing  additional  assumptions 
on  the  aa*8  (defined  in  assumption  II).  Note  also  that  assumption  (3.21)  is  always 
assumed  in  the  classical  superconvergence  theories. 

Remark  3.9.  Under  assumption  (3.21)  Theorem  1  states  that:  A  point  m  in  the. 

element  f  is  asymptotically  ri%-superconvergence  point  for  F( v)  =  —  if  and  only 

0xi 

if 


(3.22.)  e(,i;F;I;Q,x’,l,f)  <  jjL 

where 

(3.224)  e(i;F;I;<?,«',l,?)  :=  ,  *{*)■.=  |*|, 

Remark  3.10.  The  term  y%-snperconvergence  is  not  precisely  descriptive.  We 
have  used  it  because  we  wanted  that  T)  %-  super  convergence  points  for  ij  =  0  be 
superconvergence  points  in  the  common  terminology. 

Theorem  1  shows  that  the  T;%-super convergence  can  be  analyzed  by  a  local 
analysis  in  the  neighborhood  of  the  elements  under  consideration.  By  this  we 
mean  that  the  rj(%-8uper convergence  regions  of  the  finite  dement  solution  are  very 
dose  to  the  ones  determined  from  the  local  analysis  provided  that  the  mesh  is 
sufficiently  fine. 

4  The  methodology  for  determining  the  rj%- 
superconvergence  quantities 

In  this  paper  superconvergence  is  treated  as  local  behavior  and  is  based  on 
assumptions  about  the  local  behavior  of  the  solution  in  the  interior  of  the  domain. 
We  will  consider  the  class  of  solutions  which  are  locally  smooth  in  5(»°,  E),  namdy, 

(4.1)  :=  {«  6  ff!(t*)|  0<  |a|<p  +  2} 
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where  S(m°,H )  denotes  an  interior  subdomain  of  interest  in  which  the  mesh  is 
locally  periodic  as  described  above  (the  subdomain  must  be  a  finite  distance  away 
from  the  boundary  and  points  of  roughness  of  the  source  term;  see  Fig.  2a).  In 
many  instances  we  are  only  interested  in  the  subclass  of  solutions  in  UG  which  are 
“harmonic”  (we  say  that  u  is  “harmonic”  if  it  satisfies  the  homogeneous  differential 
equation),  namely, 


(4.2)  U'B "  :=  {u  €  UG  |  £( u)  =  0  in  fl} 

We  may  also  assume  that  the  functions  are  “harmonic”  in  a  subdomain  which  is 
slightly  bigger  than  S(»°,H)  and  which  includes  5(«°,ff)  in  its  interior. 

In  the  previous  Section  we  outlined  Theorem  1  which  states  that  we  can  obtain 
the  asymptotic  values  of  the  error  i/)h  for  any  smooth  solution  u  in  the  interior  of 
a  periodic  mesh-subdomain  by  solving  the  periodic  boundary-value  problem  (3.7), 
using  the  master-mesh  T  over  the  master-cell  e,  with  data  obtained  from  the  local 
(p  +  1)- degree  Taylor- series  expansion  of  the  exact  solution  (which  approximates 
well  the  function  u  in  5(a°,  H)).  Based  on  this  result  we  will  construe*  ■»  umerical 
procedure  to  determine  the  q%-Buper convergence  quantities  for  a  give-  class  of 
smooth  solutions  by  employing  the  corresponding  class  of  (p-f-l)-degree  monomials. 

a.  The  classes  of  (p  -f  1  )-degree  monomial  solutions. 

Let  us  assume  that  for  a  given  locally  periodic  grid  with  corresponding  periodic 
master-mesh  T,  given  material  orthotropy  and  given  class  of  smooth  solutions  U 
we  consider 

(4.3)  Q  :=  Iq  |  <?(*a, x 3)  =  £  ahQh(zuxi),  x\  z? 1_/} 

the  class  of  (p+  l)-degree  monomials  which  occur  in  all  (pH-  l)-degree  Taylor- series 
expansions  of  functions  from  14.  Here  Qk,  k  =  1, . . .  ,nd  denotes  a  set  of  linearly 
independent  monomials  which  span  Q.  For  example  let  us  assume  that  11  is  the 
class  of  smooth  solutions  UG  given  in  (4.1);  in  this  case  we  may  choose 

(4.4)  Qfc(xj,  s3)  :=  ,  l<fc<nd  =  p  +  2 

and  we  obtain  the  class  of  all  (p  +  l)-degree  monomials  QG. 

In  many  practical  applications  we  are  interested  only  in  the  class  of  “harmonic” 
solutions  14“ ^  given  in  (4.2).  Then  Q  is  the  two-dimensional  linear  space  of  “har¬ 
monic”  monomials  of  degree  (p  -f  1)  and  will  be  denoted  by  Q“H”.  For  example  in 
the  case  of  Laplace’s  equation  (K  =  /)  we  have 
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(4.5.)  =  . 


(4.56)  Qf  (*jf  *2)  =  TUi^1) ,  Qf(*t ,  *j)  =  Jm(*H'1) ,  z  =  z1+usa  . 

Note  that  in  this  case  we  do  not  use  the  quotes  since  the  functions  are  harmonic 
in  the  classical  sense  (i.e.  they  satisfy  Laplace’s  equation). 

6.  Freezing  the  periodicity  for  general  meshes. 

In  the  previous  Section  the  theoretical  results  for  the  asymptotic  values  of  the 
error  in  the  interior  of  periodic  mesh-subdomains  in  locally  periodic  grids  were 
outlined.  In  particular,  it  was  stated  that  the  asymptotic  error-function  yf>h  in 
any  dement  r  in  the  6-cell  in  the  interior  of  any  periodic  mesh-subdomain  can  be 
obtained  by  solving  a  periodic  boundary-value  problem  (3.7),  using  the  master- 
mesh  T  over  the  master-cell  e,  with  data  obtained  from  the  local  (p  +  1)- degree 
Taylor-series  expansion  of  the  exaat  solution.  In  order  to  apply  the  results  of  the 
theoretical  study  to  the  practical  meshes,  for  which  the  mesh  is  not  locally  periodic, 
the  following  technique  of  freezing  the  periodicity  will  be  employed: 

1.  Let  X  be  a  vertex  of  the  mesh.  We  will  be  interested  in  the  q%-superconver- 
gence  regions  in  the  dements  in  the  mesh-cell  (see  for  example  Figs.  3a, 
41)  in  a  neighborhood  of  dements  around  the  vertex  X.  Define  an  s-layered 
mesh-cell  of  dements  surrounding  the  mesh-cdl  (see  Figs.  3a  and  41 
where  <*£  is  shown  with  a  thick  black  perimeter;  this  cell  is  enclosed  in  the 
general  mesh  given  in  Fig.  1)  by 

(4.6)  uj)  :=  (J  ,  wj J:=wxDr;  u>x  :=  (J  t'; 

X€JV(t*)  X€W(t*) 


where  N(r')  denotes  the  set  of  the  three  vertices  of  dement  r1  and  u>x  is  a 
mesh-cdl  of  dements  connected  to  vertex  X  which  indudes  the  dement  r. 

2.  Complete  the  mesh-cell  ulj  (shown  in  Fig.  3b)  to  a  periodic-grid  over  a 
slightly  larger  square  periodic-cell  which  encloses  u/J,  as  shown  in  Fig.  3c. 
The  periodic-cell  is  then  scaled  and  translated  to  the  unit  master-cell  c. 
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3.  Assume  that  the  mesh  in  the  neighborhood  of  element  r  is  made  from  the  pe¬ 
riodic  repetition  of  h-cells  obtained  from  the  master-cell,  as  shown  in  Fig.  2b 
(where  the  master-cell  is  given  in  Fig.  2c),  and  let  h  tend  to  zero.  Then 
the  theoretical-setting  of  the  previous  Section  applies  and  the  asymptotic 
error-function  i>h  in  the  element  r  can  be  obtained  from  the  interpolation- 
error  and  the  solution  of  the  periodic  boundary-value  problem  (3.7)  over  the 
master-cell  c. 

In  the  numerical  examples  it  will  be  demonstrated  that  the  results  of  tj%- 
superconvergence  obtained  from  the  above  methodology  hold,  for  all  practical  pur¬ 
poses,  for  the  complex  grids  used  in  engineering  computations  (provided  that  the 
pollution-error  is  controlled,  the  approximation  is  in  the  asymptotic  range  and  that 
sufficient  number  of  mesh-layers  are  included  in  the  mesh-cell  a£). 

c.  Determination  of  the  ij%-superconvergence  quantities. 

Based  on  the  results  of  the  theoretical  study,  outlined  in  Section  3,  the  asymp¬ 
totic  locations  of  the  77%-super  convergence  quantities  in  the  element  r  for  a  solution 
«  (resp.  class  of  solutions  U)  can  be  obtained  from  the  corresponding  quantities, 
defined  on  the  master-cell  c,  for  the  local  (p  +  l)-degree  Taylor  expansion  Q  of  u 
(resp.  the  class  Q  of  the  local  (p  +  1)- degree  Taylor  series  expansions  of  functions 
it  6  U)y  namely: 


(4.7a)  Urn  K-,T,Th)  =  Cfi^Q.K^T) 

(4.76)  Um  BJ* ,(«.;  K ; t, T„)  =  B#„,(Q;  K ; f , f ) 

(4.7c)  Um  A^,(W;Jir:r,rJ  =  ^(C; *;*,?) 

(4.7<0  Um  nfM(U-,K-,r,Tk)  =  ti%u)(Q;K-,T,T) 


where  the  superconvergence  quantities  over  f  €  T  are  defined  as  the  correspond¬ 
ing  quantities  in  (2.2)  by  employing  the  relative  error  function  given  by  (3.22b). 
The  above  equalities  are  understood  through  the  transformation  which  maps  the 
master-cell  into  the  /i- cells. 

Remark  J^.l.  The  above  limits  hold  for  the  locally  periodic  meshes  under  the 
assumptions  of  the  theoretical  analysis  of  Section  3.  Hence  for  the  general  grids 
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the  limit  should  be  understood  for  the  mesh  which  is  constructed  by  freezing  the 
periodicity  (see  also  Figs.  2  and  3). 

The  asymptotic  i/%- contours  for  a  given  solution  u  can  be  obtained  by  con¬ 
touring  the  function  defined  in  (3.8),  corresponding  to  the  local  Taylor- series 
expansion  Q  of  «.  The  superconvergence  points  I  for  a  given  class  of  solutions  U 
should  satisfy 

(4.8)  F(^ )(!)  =0  for  1  <  *  <  nd 


Here  ^  :=  pt  —  zf ,  where  sf  is  the  solution  of  (3.7)  for  p{  —  Q{  —  (Qi)?iT  where  Q{ 
is  the  i-th  basis  monomial  of  the  nd-dimensional  monomial  space  Q  corresponding 
to  the  class  U ,  as  discussed  in  (4.3),  (4.4).  When  (4.8)  holds,  the  zero-contours  of 
Ffa)  intersect  at  i  for  1  <  i  <  nd. 

Remark  4.2.  For  the  class  of  general  solutions  Z/°,  nd  =  p  +  2,  and  as  p  increases 
(for  p  >  3)  the  likelihood  of  zero-contours  for  each  one  of  the  ^ ’s  crossing  simul¬ 
taneously  at  one  point  is  very  small.  Thus  a  few  or  no  superconvergence  points 
exist  for  the  class  U°  and  p  >  3.  On  the  other  hand  for  the  class  of  “harmonic” 
solutions  U"H' ,  nd  =  2  independent  of  p,  and  many  super  convergence  exist  and 
their  number  increases  with  p.  For  details  see  [28]. 

The  asymptotic  17%- superconvergence  bands  for  a  solution  «  can  be  determined 
from  the  function  F(V>)  by  using  piecewise  linear  interpolation  of  F(V»)  on  a  suffi¬ 
ciently  fine  uniform  mesh  obtained  by  subdividing  the  triangle  r. 

The  asymptotic  17%-superconvergence  regions  for  a  class  of  solutions  U  can  be 
determined  by  using  numerical  optimization.  In  particular,  consider  the  uniform 
subdivision  of  the  element  r  into  subtriangles  with  vertices  at  the  set  of  points 
5  :=  Define  the  relative-error  function  at  the  point  £k 

lEa^XOI 

(4.9)  efoo&ki  F)  K;  c,  {*?}&,  l,f)  :=  max - - 100 , 

tSl 

Then  the  function  F;  K\  G,  1>  ?)  can  be  defined  for  any  point 

m  €  f  by  using  linear  interpolation  in  the  subtri angles.  The  asymptotic  rj%- 
superconvergence  regions  in  the  element  r  can  be  approximated  using  the  level-sets 
of  the  functions  6f(w)(«;  F;  K\  Q,  {zfJSfi,  l,f)  i.e. 

(4.10)  K-,f,T) »  {*  €  f  |  FiiTie,{«fRf1, l.f)  <  v%] 
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We  will  call  thii  approach  the  dirtct  approach.  It  it  also  possible  to  use  a  timpUfied 
approach  which  avoids  the  use  of  numerical  optimisation  at  every  point.  First 
observe  that 


(4.11) 

i£aiF’(^i)(i)| 

max - — _ 

.*“*  |  £  a,-  m  )(£,)  | 

tael 


max 

•< 


(E«?)* 


|£ojF(*i)(i)|'V 

is] 


V  J?"L  I  E  ^  )(^)l  (  £  a?)* 

tsl  tsl 


<  max 


(E  «?)» 

tsl 


l<  jnax  1 2  (*.)({,)! 

isl 


-g(w,x*)),)‘< 


1 2«<wi)(^)l] 

•  J“WS  tael  '  I 


nun 

a. 


(E  «?)* 

i=l 


(|(mx*))1) 


Hence  we  can  define  the  function 


(4.12)  JTsQ. {O&.l.f)  :=  ^  ^g(*(*X*))* . 


where 


(4.13) 


Za:= 


mm 


j=l,...,wp  i=l _ _ _ 

(E«?)* 

4  1 


The  quantity  ZE  can  be  computed  using  numerical  optimisation.  Let 


(4.14)  :=  {•  6  r  I  §!<.,(*; F-.K-.Q, K}&, l,f)  <  ^  } 

denote  the  approximate  regions  of  q%~super  convergence  for  the  class  of  solutions 
Q  obtained  by  the  simplified  approach. 
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The  common  q%-superconvcrgcnce  region  in  an  dement  r  for  a  given  class  of 
meshes  and  material!  can  be  determined  by  using  either  the  direct  or  the  simplified 
approach. 

Remark  4-8.  Note  that  we  have 
(4.18)  l%u)(Q-,K-,t,i)  c 

and  we  get  a  conservative  estimate  for  the  17%-superconvergence  regions  when  the 
amplified  approach  is  used.  For  the  class  of  harmonic  functions  U’B’  we  observed 
that  the  regions  obtained  by  the  two  approaches  are  very  close.  However  for  the 
class  of  general  solutions  UG  the  simplified  approach  is  too  pessimistic  and  the 
more  expensive  direct  approach  must  be  employed. 

Remark  4.4.  The  functions  defined  in  (4.9),  (4.12)  depend  on  the  set  of  points  E. 
To  ensure  good  accuracy  in  the  approximation  of  the  q%-  super  convergence  regions 
a  sufficient  number  of  points  must  be  employed. 

5  Numerical  studies  of  superconvergence 

We  will  now  present  a  model  study  of  the  17  %-  superconvergence  properties  of 
finite  dement  solutions  in  the  interior  of  any  grid,  for  any  smooth  solution  and 
material-orthotropy.  We  will  address  the  following  questions: 

1.  Given  an  dement  r  in  the  interior  of  a  periodic  subdomain  of  a  locally  pe¬ 
riodic  mesh  T*.  where  are  the  points  of  superconvergence  and  what  is  their 
sensitivity  with  respect  to: 

(i)  The  class  of  solutions  (harmonic  or  general); 

(ii)  The  material  orthotropy; 

(iii)  The  geometry  of  the  mesh  and  the  number  of  mesh-layer  in  the  mesh- 
cdlswj? 

2.  Given  q%  (between  0%-100%)  where  are  the  regions  in  which  the  derivatives 
of  the  error,  or  the  components  of  the  flux  are  within  17%  of  the  maximum 
error  in  the  respective  quantities  in  the  dement  r?  In  particular: 

a.  For  the  periodic  meshes  shown  in  Fig.  4  where  are  the  q  %-  super  conver¬ 
gence  regions  for  the  ^-derivative  and  how  do  these  regions  vary  with 
the  degree  p  of  the  dements: 

(i)  For  the  class  of  harmonic  solutions  UB\ 
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(ii)  For  the  class  of  general  solutions  Ual 

b.  For  periodic  meshes  what  is  the  effect  of  the  mesh- distortion,  the  material- 
orthotropy  and  the  surrounding  mesh  on  the  location  of  the  17%-supercon- 
vergence  regions  in  a  typical  element? 

c.  For  a  general  mesh,  like  the  ones  used  in  practical  engineering  compu¬ 
tations,  how  different  are  the  regions  which  are  determined  by  freezing 
the  periodicity  from  the  regions  determined  from  the  actual  mesh? 

We  will  answer  the  above  mentioned  questions  by  employing  the  computer- 
based  approach  outlined  in  Section  4. 

5.1  Robustness  of  the  superconvergence  points  with  re¬ 
spect  to  the  class  of  solutions  and  distortions  of  the 
mesh  for  the  Laplace  equation  (K  =  I) 

We  considered  the  four  periodic  patterns  shown  in  Fig.  4.  In  [28]  we  determined 
the  points  of  superconvergence  for  the  periodic  meshes  shown  in  Fig.  4  for  the 
dasses  UB  and  U°.  The  points  are  listed  in  Tables  1  for  1  <  p  <  5  for  the  class 
of  harmonic  solutions.  Table  2  indicates  which  points  in  each  pattern  are  also 
superconvergence  points  for  the  dass  of  general  solutions.  In  order  to  investigate 
the  robustness  of  these  points  with  respect  to  the  dass  of  solutions  we  computed 
(Table  3)  the  maximum  q%- error  in  the  xr derivative  of  the  solution  for  the  dass  of 
general  solutions  UG  at  the  superconvergence  points  for  the  dass  UH  for  quadratic 
dements  i.e.  the  value 


(5.1)  i7*(i)  :=  ,  *  e 

From  the  numerical  results  we  observe  that  the  rdative  error  at  the  supercon¬ 
vergence  points  for  the  dass  of  harmonic  solutions  the  may  be  large  for  solutions 
which  are  not  harmonic. 

5.2  Determination  of  the  17%-superconvergence  regions  for 

the  periodic  meshes  (K  =  I) 

From  the  previous  example  it  is  dear  that  the  superconvergence  points  are 
very  sensitive  to  dass  of  solutions  and  the  geometry  of  the  mesh.  Moreover,  the 
points  may  not  exist  at  all;  for  example  there  are  no  superconvergence  points  for 
the  -derivative  of  the  solution  in  the  Union-Jack  pattern  for  the  dass  of  general 
solutions.  We  will  now  report  the  q%-super  convergence  regions  for  the  periodic 
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meshes  shown  in  Fig.  4.  We  will  show  that,  unlike  the  superconvergence  points, 
these  regions  always  exist  (in  any  element  in  any  pattern)  for  sufficiently  high  17%. 

We  determined  the  77%-super  convergence  regions  for  the  classes  of  harmonic 
and  general  solutions,  for  the  elements  f  in  the  periodic-cell  which  consists  of:  two 
dements  for  the  Regular  pattern;  four  elements  for  the  Chevron  and  Criss-Cross 
patterns;  eight  dements  for  the  Union- Jack  pattern.  The  levels  of  17%  in  the  Figures 
were  chosen  to  be  10%,  20%,  30%.  In  Figs.  5-8  the  regions  and 

the  contours  of  zero-error  in  the  x1-derivative  of  the  solution  for  tke  basis  harmonic- 
monomials  are  given  for  p  =  1,  2,  3,  4,  5.  It  should  be  noted  that  for  the  dass  of 
harmonic  solutions  the  results  were  obtained  using  the  “simplified  approach”  (these 
results  are  very  dose  to  corresponding  results  obtained  using  the  direct  approach). 
The  regions  iV1]*  (Q?;I;f,f)  for  p  =  1,  2,  3,  are  given  in  Figs.  9-12.  It  should  be 

noted  that  for  the  dass  of  general  solutions  the  direct  approach  must  be  used  in 
the  determination  of  the  17%-superconvergence  regions,  as  the  simplified  approach 
is  too  pessimistic  (it  results  in  regions  which  are  much  smaller  than  the  actual 
ones). 

We  observe  that: 

1.  For  the  dass  of  harmonic  solutions,  the  17%-superconvergence  regions  for  the 
error  in  the  -derivative  grow  with  the  degree  p  of  the  dements. 

2.  For  the  dass  of  general  solutions,  the  17%-superconvergence  regions  for  the 
Xj -derivative  are  small  (except  for  the  Regular  pattern  for  p  =  1,  2)  and  they 
diminish  in  size  as  p  grows,  i.e.  for  higher  p’s  (p  >  3)  the  regions  may  not 
exist  at  all  for  small  17’s,  for  example  for  17  <  30. 


5.3  Robustness  of  T7%-superconvergence  regions:  Effect  of 
mesh-distortion  (K  =  I) 

To  study  the  robustness  of  the  superconvergence  points  with  respect  to  dis¬ 
tortions  of  the  mesh  we  considered  the  Regular  pattern  with  the  central  node 
displaced  (see  Fig.  13)  and  computed  the  maximum  17%-error  in  the  x1- derivative 
of  the  solution  at  the  points  with  the  same  master-dement  coordinates  as  the 
superconvegence  points  in  the  undistorted  mesh  (given  in  Table  la)  i.e. 


(5.2) 

«*(*')  • 


where  f  =  T{f)  and  T  :  T  — ►  T'  denotes  the  transformation  from  the  undistorted 
master-mesh  T  to  the  distorted  master-mesh  T'  and  denotes  the 
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set  of  super  convergence  points  in  the  undistorted  pattern.  The  value  of  f}%  was 
computed  for  both  the  class  of  harmonic  and  class  of  general  solutions  (Table  4), 
for  quadratic  elements. 

In  order  to  study  the  robustness  of  the  T;%-Buper convergence  regions  with  re¬ 
spect  to  distortions  of  the  mesh,  the  regions  'fcj*  (QB\  J;  f ,  T)  were  determined  as 

a  function  of  the  mesh-distortion  for  periodic  grids  of  quadratic  elements  in  the 
Regular  and  Union-Jack  patterns.  In  both  patterns  the  central  node  was  moved  to 
(0.4,  0.4),  (0.6,  0.6),  (0.8,  0.9)  from  its  position  (0.5,  0.5)  in  the  undistorted  pat¬ 
tern,  as  shown  in  Fig.  13.  The  regions  7T&  (Qb;I;t,T)  (if  =  10, 20, 30)  are  given 

in  Figs.  14,  15  (Fig.  14  for  the  Regular  pattern,  and  Fig.  15  for  the  Union-Jack 
pattern). 

We  make  the  following  observations: 

1.  For  both  the  classes  of  harmonic  and  general  solutions,  the  17%-error  at  the 
points  with  master-coordinates  corresponding  to  the  superconvergence  points 
in  the  undistorted  pattern  can  vary  from  2.5%  to  100%  when  the  pattern  is 
distorted. 

2.  For  small  distortions  of  the  pattern,  i.e.  when  the  central  node  is  located  at 
(0.4, 0.4)  or  at  (0.6, 0.6)  the  regions  for  17%- error  (if  =  10, 20, 30)  do  not  differ 
much  from  the  regions  in  the  undistorted  pattern  and  there  are  big  common 
regions  for  all  the  distorted  meshes  and  the  values  of  if%  considered. 

3.  For  large  distortions  of  the  pattern,  i.e.,  when  the  central  node  is  located  at 
(0.8, 0.9)  the  common  regions  for  all  the  distorted  meshes  exist  only  for  large 
values  of  q%  (if  >  30). 

5.4  Robustness  of  ?7%-superconvergence  regions:  Effect  of 
grid- material  orientation  (K  ^  J) 

Given  a  mesh  and  an  orthotropic  material,  we  would  like  to  find  the  i)%- 

superconvergence  region  ^jF>*u)(G"F";  M;  f,  T)  for  F( u)  =  or  q^u)  which  is 

common  to  a  given  class  of  orthotropic  materials  A i.  We  considered  the  Union- Jack 

pattern  and  quadratic  elements  and  let  q%  =  25%.  We  assumed  that  =  2  or 

4  and  we  let  angle  9  of  orientation  of  the  principal  plane  of  orthotropy  with  respect 
to  the  mesh  vary  horn  0°  to  90°.  We  use  the  notation 

(5.3)  Mr  :=  {*  =  K(9)\  — ^  =  r  and  9  €  [0,90°]} 

l  A  sain  * 
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We  determined  the  common  regions  £**)(£"jr>;  M3;  f,  f)  (shown  in  Fig.  16)  and 

(ghown  in  Fig.  17).  The  common  25%-super  convergence  re¬ 
gions  were  approximated  by  finding  the  intersections  of  the  25%-super  convergence 
regions  for  9  =  0°,  15°,  30°,  45°,  60°,  75°,  90°  as  shown  in  Figs.  16,  17,  18.  In 
Fig.  18  we  give  the  common  region  ;  Af 4;  f ,  T).  It  can  be  seen  that 

jr  1 

for  T7~—  —  2  the  common  regions  of  25%  error,  for  the  Zj-derivative  or  the  xx- 

Amm 

components  of  the  flux  for  all  6  €  [0, 90°],  are  of  significant  size  with  respect  to  the 
size  of  element.  However,  for  -  m*x  =  4  the  common  regions  of  25%  error,  for  the 

Amtii 

-derivative  for  all  6  €  [0,90°],  are  small.  It  is  clear  that  the  25%-common  region 
will  dissapear  for  higher  values  of  the  ratio  — ' 

Amin 


5.5  Robustness  of  the  ^%-superconvergence  regions: 

Effect  of  the  surrounding  mesh 

Here  we  shall  demonstrate  that  the  »j%-Buper  convergence  regions  and  contours 
in  an  element  r  are  affected  only  by  one  or  two  surrounding  layers  of  elements. 
We  considered  five  3  x  3-cell  periodic  meshes  Tx  —  Tt ,  as  shown  in  Fig.  19.  We 
determined  the  regions  )  i  —  1,.  ..,5,  where  ij  =  10,  20,  30  and 

for  the  elements  of  degree  p  =  2,  3,  4.  The  f/%-superconvergence  regions  and 
the  0%  contours  are  shown  in  Figs.  20-22.  Fig.  23  displays  the  common  region 
J;f,X),  where  T  =  {Tt  }?=1,  (shaded  gray)  for  the  five  3  x  3-cell  meshes 

for  p  =  2, 3.  Next,  we  added  one  layer  of  elements  around  the  five  3  x  3-cell  meshes, 
to  get  the  five  5  x  5-cell  meshes  17,  i  =  1, ...  ,5,  as  shown  in  Fig.  24.  Figs.  25-27 
display  the  regions  (Qa]I;f,T-),  i  =  1,...,5,  for  p  =  2,  3,  4.  In  Fig.  28 

we  show  the  common  region  7^*jff(Qg;I;  where  V  —  for  the  five 

5  x  5-cell  meshes  for  p—  2,  3. 

From  the  results  we  observe  that: 

1.  The  regions  *E^  (QB\  I ;  f,  T{ ),  are  essentially  identical  with  the  correspond¬ 
ing  regions  rtg.  ( QB ;  J;  f ,  T-).  However,  the  zero-contours  for  the  Zj- derivative 

of  the  error  for  the  harmonic  monomials  are  different  for  the  corresponding 
grids  Tt ,  T/,  *  =  1,. . .  ,5. 

2.  The  common  regions  7t*^(Qg;I;f,T)  and  7lij£(Q,a]I;T,T')  are,  for  all 
practical  purposes,  identical. 
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To  see  whether  the  mesh  beyond  two  mesh-layers  of  elements  affects  the  com¬ 
mon  ?;%-8uper convergent,-  region  in  element  f  and  the  zero-contours  for  the  zx- 
derivative  of  the  error  for  the  basis  monomials,  we  surrounded  the  5  x  5-cell  mesh  of 
Fig.  24e  by  5  more  layers  of  uniform  elements  to  get  the  mesh  Tg'  shown  in  Fig.  29. 
From  Fig.  30  we  can  see  that  the  t;%- super  convergence  regions  TVet  (QH\I\r,T") 

for  the  element  f  are  identical  to  the  corresponding  contours  and  q%-superconvergence 
regions  for  the  5  x  5-cell  mesh  Tg.  We  also  surrounded  the  5  x  5-cell  mesh  of  Fig.  24e 
by  5  layers  of  distorted  elements  to  obtain  the  mesh  T5'",  as  shown  in  Fig.  31.  From 
Fig.  32  we  can  observe  that  the  regions  for  the  element  f  are 

practically  identical  to  those  for  the  5  x  5-mesh  Tg. 

From  the  above  examples  we  can  conclude  that:  The  i}%-superconvergence  re¬ 
gions  in  element  f  are  not  affected,  for  all  practical  purposes,  by  the  mesh  outside 
the  immediate  £  layers  of  elements  surrounding  the  element. 

5.6  77%-superconvergence  for  general  meshes 

In  order  to  demonstrate  that  the  method  of  freezing  the  periodicity  gives 
accurate  prediction  of  the  q%-  super  convergence  for  general  meshes  we  considered 
the  following  examples: 

1.  An  example  which  shows  how  the  conclusions  from  the  periodic  grid  translate  to 
grids  in  finite  domains. 

We  considered  a  5  x  5  mesh  of  Union-Jack  cells  surrounded  by  an  arbitrary 
quasiuniform  mesh  generated  by  a  remeshing  algorithm.  For  the  mesh  and  the 
mesh- cells  shown  in  Figs.  33a-33c,  we  determined  the  t/%- superconvergence  re¬ 
gions  for  the  xx -derivative  of  the  error  for  cubic  elements  and  the  class  of  quartic 
harmonic  polynomial  solutions.  In  Figs.  34a-34c  we  give  the  17%-superconvergence 
regions  for  the  xx -derivative  for  the  mesh-cells  shown  (shaded  gray)  in  Figs.  33a- 
33c,  respectively.  These  regions  were  determined  from  the  approximate  solutions 
of  Dirichlet  boundary-value  problems  over  the  finite-domain  using  data  obtained 
from  the  basis  harmonic-monomials.  Comparing  Figs.  34a-34c  (for  the  Union-Jack 
cells  in  the  mesh  which  covers  the  finite-domain)  with  Fig.  7c  which  gives  the  17%- 
superconvergence  regions  for  the  periodic  Union-Jack  cell,  we  can  observe  only  a 
slight-difference  between  the  corresponding  regions  shown  in  Fig.  7c  and  Fig.  34c. 
This  example  indicates  that  the  results  on  the  q%-superconvergence  regions  ob¬ 
tained  from  the  periodic  cell  hold,  for  all  practical  purposes,  even  for  cells  very 
close  to  the  boundary  of  the  periodic  submesh. 

£.  An  example  which  demonstrates  that  the  results  obtained  using  the  method 
of  freezing  the  periodicity  hold  for  arbitrary  meshes  and  solutions  which  are  not 
polynomials. 
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We  considered  the  mesh  shown  in  Fig.  1  in  the  square  domain  Cl  =  (0,  l)3  and 
identified  several  elements  in  small  mesh-cells  enclosed  in  the  interior  of  the  mesh. 
The  mesh-cells  uj}  are  shown  in  Figs.  35a-35c;  in  these  Figures  the  mesh-cells  are 
shaded  gray.  For  each  mesh-cell  we  isolated  a  corresponding  submesh  u/3  of  the 
original  grid  which  includes  the  mesh-cell  and  three  mesh-layers  around  it;  the 
submeshes  w,  are  shown  within  the  thick-line  perigram  in  Figs.  35a-35c.  Each  one 
of  the  submeshes  was  scaled  and  completed  to  a  square  periodic-cell  c,  as  shown 
in  Figs.  36a-36c.  We  considered  the  class  of  the  following  four  harmonic  solutions: 


(5.4a)  Ui(rit0i )  -  r f*  sin(a,^) ,  *  =  1, ...  ,4 , 


(5.46)  r?  =  (*!  -  *i)3  +  (*j  -  z*j)J ,  6i  =  tan"1  (**--*?) 
The  solutions  are  specified  by  the  following  parameters: 


(5.5) 


(xiixa)=(2.,0.)  ,  <*1  =  ^,  (*i>*a)  =  (1-5, 0.5)  ,  a2  =  |, 


[  (xi»xa)  —(0.5,2.) ,  as  —  n  >  (xx>xa)  -  (  3*»0.)  ,  <*4  —  7  • 


We  determined  the  common  region  j&jpjU;  I ; r,  Th)  from  the  actual  mesh  and 
'R}q?{U\  I ;  f ,  T)  from  the  periodic  meshes  T  obtained  by  freezing  the  periodicity  lo¬ 
cally,  U  =  {iii}J=1  for  linear  and  quadratic  elements.  The  regions  ’Ea|*(W;  J;t,T) 

were  obtained  by  solving  Neumann  boundary-value  problems,  using  data  from 
the  exact  solutions,  over  the  entire  domain  Cl  and  using  the  mesh  Th  of  Fig.  1. 
The  regions  7£"|?(W;J;t,T)  were  obtained  by  solving  periodic  boundary-value 

problems,  using  data  obtained  from  the  local  Taylor-series  expansion,  over  the 
periodic-cell  and  the  periodic  meshes  of  Fig.  36a-36c.  In  Figs.  37a-37c  we  give 
the  region  7£3|?(W; I;r,Th)  for  linear  elements  while  Figs.  38a-38c  we  show  the 


corresponding  regions  'R?^(U\  I]f,T)  obtained  from  the  solutions  of  the  periodic 

boundary-value  problems.  In  Figs.  39a-39c  (resp.  Figs.  40a-40b)  we  give  the  region 
*’«(«; l-,r,Th)  (resp.  J;f,T))  for  quadratic  elments. 
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We  observe  that  the  regions  jt*g(U’, IjrtTh)  obtained  by  using  the  solutions 

of  the  Neumann  boundary- value  problems  over  the  entire  grid  practically  coincide 
with  ,  the  regions  t,  T)  determined  from  the  periodic  boundary- value 

problems  over  the  periodic  meshes.  Hence  the  method  of  freezing  the  periodicity 
gives  results  which  hold  for  the  practical  grids. 

5.7  Robustness  of  the  q%-superconvergence  points  in  prac¬ 
tical  computations. 

Here  we  give  an  example  which  shows  that  T} %-  superconvergence  points  are 
robust  and  can  be  used  to  sample  the  solution  in  the  interior  of  any  grid  used  in 
practical  computations.  For  the  mesh-cell  shown  shaded  gray  in  Fig.  41,  we 
determined  the  (fig;  for  i)  =  10,20,30,  for  p  =  1,  2,  3.  The  i}%- 

regions  are  shown  in  Fig.  42a,  42b,  42c  for  p  =  1,2,  3,  respectively.  We  extracted 
some  random  sampling  points  within  the  10%-regions  for  the  a;  j- derivative  of  the 
solution  for  the  class  of  harmonic  solutions  (as  shown  in  Fig.  43)  and  we  solved  a 
Neumann  boundary  value  problem  over  the  mesh  shown  in  Fig.  1,  with  the  exact 
solution  ttj  and  domain  fl*  given  as: 

(  )-  1  1  1  1  1 

1.5a  —  (z  +  Zo)a  1.5*  +  (z  +  ZoY  1.5a  -  1.0  +  1.5a  +  1.0 

where  Zq  =  0.5(1  + 1) ,  fl1  =  (— 0.5, 0.5)3,  and 


(5.65)  «,(*!,*,)  =  fte((z  -  *„)*) 

where  z0  =  0.5  +  2.0*,  ft2  =  (0,  l)3. 

For  the  two  solutions  given  above  we  determined  the  relative  error  in  the  z1- 
derivative  of  the  solution  at  the  sampling  points  for  p  =  1,  2,  as  shown  in  Table 
5a,  5b.  From  the  results  we  observe  that: 

(1)  As  p  increases,  the  rj%  regions  for  the  zx -derivative  of  the  solution  increase 
in  size  relative  to  the  elements. 

(2)  From  the  results  shown  in  Table  5a,  5b,  we  can  see  that  the  relative  error  at 
the  sampling  points  is  less  than  15%  for  the  two  solutions  considered. 

Thus,  this  example  shows  that  from  the  17%-regions  obtained  from  the  computer- 
based  approach  of  this  paper  we  can  extract  robust  sampling  points ,  which  can  be 
used  for  any  solution  in  the  given  class. 
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6  Summary  of  conclusions 

1.  The  new  concept  of  17 %-super convergence  was  introduced  and  the  17%- super¬ 
convergence  regions  for  various  flux-quantities  and  for  various  classes  of  so¬ 
lutions  and  materials  were  defined. 

2.  The  77  %-super  convergence  regions  can  be  found  for  any  mesh  using  a  com¬ 
puter  based-approach. 

3.  The  q%-8uper convergence  regions  for  the  Xj- derivative  were  determined  for 
four  periodic  patterns  and  it  was  observed  that: 

a.  The  regions  always  exist  for  sufficiently  high  17. 

b.  For  fixed  value  of  17  the  regions 

(i)  always  exist  for  the  class  of  harmonic  solutions  and  grow  in  size 
relative  to  the  element  as  the  degree  p  of  the  elements  is  increased; 

(ii)  may  not  exist  for  the  class  of  general  solutions  and  shrink  in  size  as 
the  degree  p  of  the  elements  is  increased. 

These  observations  can  be  expected  to  hold  for  ^  riodic  grids  made  by  the 
repetition  of  any  periodic  pattern. 

4.  The  17%-supercon vergen ce  regions  can,  for  all  practical  purposes,  be  deter¬ 
mined  for  any  general  finite  element  grid,  like  for  example,  the  grids  employed 
in  practical  engineering  computations  by  using  the  technique  of  freezing  the 
periodicity.  Thus  the  new  concept  of  77%- superconvergence  can  be  employed 
to  identify  optimal  sampling  points  for  the  flux-quantities  in  the  interior  of 
any  practical  grid. 

5.  We  defined  the  common  77  %-super  convergence  regions  for  entire  classes  of 
grids  and  materials.  These  regions  can  be  determined  numerically  and  can 
be  used  to  identify  robust  sampling  points  for  the  flux-quantities  of  interest 
for  entire  classes  of  grids,  solutions  and  materials. 
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List  of  Figures 

Fig.  1.  Typical  example  of  a  general  finite-element  grid  used  in  practical  engi¬ 
neering  computations.  This  grid  was  generated  using  a  commercial  mesh-generator 
for  a  problem  with  a  material  interface. 

Fig.  3.  (a)  An  example  of  a  domain  with  a  locally  periodic  subdomain  in  its 

interior  which  is  covered  exactly  by  a  periodic  array  of  h- cells;  (b)  A  detailed 
depiction  of  the  mesh  Th(m°,H)  in  the  interior  of  the  periodic  mesh-subdomain 
5(sj°,  H).  This  grid  is  made  by  the  periodic  repetition  of  the  h- scaled  master  mesh 
T;  (c)  The  master  mesh  T  in  the  master-cell  c. 

Fig.  3.  Extraction  of  a  mesh  cell  w*  and  completion  to  a  periodic  mesh  T  in  the 
periodic  cell  c  (a)  The  grid  (of  Fig.  1)  with  the  mesh-cell  (with  its  perimeter 
shown  in  thick  black  line)  and  the  mesh- cell  (shaded  gray);  (b)  The  mesh- cell 

with  the  mesh-cell  <*£  (shaded  gray)  in  its  interior;  (c)  The  scaled  mesh- cell  u/* 
embedded  in  a  periodic  mesh  T. 

Fig.  4.  Periodic  meshes  of  triangles  with  the  elements  r  shaded  gray:  (a)  Regular 
pattern;  (b)  Chevron  pattern;  (c)  Union-Jack  pattern;  (d)  Criss-Cross  pattern. 


Fig.  5.  17%-superconvergence  regions  for  periodic  meshes:  The  regions 
for  the  levels  tj  =  10,  20, 30,  and  the  zero-contours 


C°^(Qf ; JT;f,T),  *  =  1*2,  for  the  Regu’ar  pattern,  (a)  p  =  1;  (b)  p  -  2;  (c)  p  =  3; 
(d)  p  —  4;  (e)  p  =  5.  The  regions  were  determined  using  the  simplified  approach. 


Fig.  6.  »7%-super convergence  regions  for  periodic  meshes:  The  regions 
fcg  (Qh;I;t,T)  for  the  levels  17  =  10,  20,  30,  and  the  zero-contours 

&£XQf  ;I;f,T),  *  =  1,2,  for  the  Chevron  pattern,  (a)  p  =  1;  (b)  p  -  2;  (c)  p  =  3; 
(d)  p  =  4;  (e)  p  =  5.  The  regions  were  determined  using  the  simplified  approach. 


Fig.  7.  17 %-super convergence  regions  for  periodic  meshes:  The  regions 
T)  for  the  levels  r\  =  10,  20,  30,  and  the  zero-contours 

Cg  i  =  1,2,  for  the  Union-Jack  pattern,  (a)  p  =  1;  (b)  p  =  2; 

(c)  p  =  3;  (d)  p  =  4;  (e)  p  =  5.  The  regions  were  determined  using  the  simplified 
approach. 


Fig.  8.  17%-Buperconvergence  regions  for  periodic  meshes:  The  regions 
(gg;J;f,[f),  for  the  levels  17  =  10,  20,  30,  and  the  zero-contours 
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*  =  1,2,  for  the  Criss-Cross  pattern,  (a)  p  =  1;  (b)  p  =  2; 

(c)  p  =  3;  (d)  p  =  4;  (e)  p  =  5.  The  regions  were  determined  using  the  simplified 
approach. 

Fig.  9.  17%-superconvergence  regions  for  periodic  meshes:  The  regions 
*3L(Q®;I;  f  ,T),  for  the  levels  t/  =  10, 20,  30,  for  the  Regular  pattern,  (a)  p  =  1; 
(b)  p  =  2;  (c)  p  =  3.  The  regions  were  determined  using  the  direct  approach. 

Fig.  10.  17%-superconvergence  regions  for  periodic  meshes:  The  regions 
fcjhi  (Q°i  J»  f ,  T),  for  the  levels  17  —  10, 20,  30,  for  the  Chevron  pattern,  (a  )p—  1; 

(b)  p  =  3.  It  should  be  noted  that  for  p  =  2  there  is  no  region  for  t/%  <  30%.  The 
regions  were  determined  using  the  direct  approach. 

Fig.  11.  17%-superconvergence  regions  for  periodic  meshes:  The  regions 

I;  f,  T),  for  the  levels  if  =  10, 20, 30,  for  the  Union-Jack  pattern,  (a)  p  = 
2;  (b)  p  =  3.  Note  that  for  p  =  1  there  is  no  region  if%  <  30%. 

Fig.  12.  17 %-super convergence  regions  for  periodic  meshes:  The  regions 
Rjk  (Qa;  I ;  f ,  T),  for  the  levels  17  =  10, 20, 30,  for  the  Criss-Cross  pattern,  (a)  p  = 
I;  O’)  P  —  2;  (c)  p  =  3.  The  regions  were  determined  using  the  direct  approach. 

Fig.  13.  Robustness  of  17%-superconvergence  regions:  Effect  of  mesh  distortion. 
The  undistorted  patterns  are  given  in  the  domain  (0,  l)x(0,  1)  with  the  central 
node  located  at  (0.5,  0.5).  (a)  Regular  pattern  with  the  central  node  at  (0.4,  0.4); 
(b)  Regular  pattern  with  the  central  node  at  (0.5,  0.5);  (c)  Regular  pattern  with 
the  central  node  at  (0.6,  0.6);  (d)  Regular  pattern  with  the  central  node  at  (0.8, 
0.9);  (e)  Union-Jack  pattern  with  the  central  node  at  (0.4,  0.4);  (f)  Union-Jack 
pattern  with  central  node  at  (0.5,  0.5);  (g)  Union-Jack  pattern  with  the  central 
node  at  (0.6, 0.6);  (h)  Union-Jack  pattern  with  the  central  node  at  (0.8,  0.9). 

Fig.  14.  Robustness  of  17%-superconvergence  regions:  Effect  of  mesh-distortion. 
Distorted  Regular  pattern,  quadratic  elements  and  cubic  harmonic  polynomial  so¬ 
lutions.  The  regions  'Rj*  {QH\  J;f,T)  when  the  central  node  is  at:  (a)  (0.4,  0.4); 

(b)  (0.5,  0.5);  (c)  (0.6, 0.6);  (d)  (0.8,  0.9).  The  regions  were  determined  using  the 
direct  approach. 

Fig.  15.  Robustness  of  17%-superconvergence  regions:  Effect  of  mesh-distortion; 
Distorted  Union-Jack  pattern,  quadratic  elements  and  cubic  harmonic  polynomial 
solutions.  The  regions  'Rj*  (QH\  I]f,T)  when  the  central  node  is  at:  (a)  (0.4, 

0.4);  (b)  (0.5, 0.5);  (c)  (0.6, 0.6);  (d)  (0.8,  0.9).  The  regions  were  determined  using 
the  direct  approach. 
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Kg.  16.  Robustness  of  17%-superamvergence  regions:  Effect  of  the  grid-material 
orientation;  Union-Jack  pattern,  quadratic  dements  and  cubic  “harmonic”  poly¬ 
nomial  solutions.  The  common  region  for  ~  =  2  and  all 

grid-material  orientations  8  €  [0°, 90°]. 

Fig.  17.  Robustness  of  t7%-super  convergence  regions:  Effect  of  the  grid-material 
orientation;  Union-Jack  pattern,  quadratic  dements  and  cubic  “harmonic”  poly¬ 
nomial  solutions.  The  common  region  for  —  2  and  all 

grid-material  orientations  8  €  [0°,90°]. 

Fig.  18.  Robustness  of  17%-superconvergence  regions:  Effect  of  the  grid-material 
orientation;  Union-Jack  pattern,  quadratic  dements  and  cubic  “harmonic”  poly¬ 
nomial  solutions.  The  common  region  'k^(QB \ M.a\t,T)  for  =  4  and  all 

grid-material  orientations  8  €  [0°,90°]. 

Fig.  10.  Robustness  of  the  17%- super  convergence  regions:  Effect  of  the  surround¬ 
ing  mesh;  The  five  3  x  3-cell  meshes  used  with  the  dement  f  subdomain  shaded 
gray,  (a)  Mesh  Tx ;  (b)  Mesh  Ta;  (c)  Mesh  Ta ;  (d)  Mesh  T4;  (e)  Mesh  T% . 

Fig.  20.  Robustness  of  the  »/%- super  convergence  regions:  Effect  of  the  sur¬ 
rounding  me 
regions 

As  —  1, . . . ,  5, 

Tj ;  (b)  Mesh  Ta ;  (c)  Mesh  Ta ;  (d)  Mesh  T4 ;  (e)  Mesh  T, . 

Fig.  21.  Robustness  of  the  17%-superconvergence  regions:  Effect  of  the  sur¬ 
rounding  mesh;  Cubic  dements,  quartic  harmonic  polynomial  solutions.  The  re¬ 
gions  £?gt(Qir;I';r,Tfc)  and  the  zero-contours  C^(Qf ;  J;f,Tfc),  t  =  1,  2,  for 

k  =  1, ...  ,5,  in  the  dement  f  for  the  five  3  x  3-cell  meshes  of  Fig.  19.  (a)  Mesh 
Tj ;  (b)  Mesh  Ta ;  (c)  Mesh  T, ;  (d)  Mesh  T4 ;  (e)  Mesh  T5 . 

Fig.  22.  Robustness  of  the  »7%-8uper convergence  regions:  Effect  of  the  sur¬ 
rounding  mesh;  Quartic  dements,  quintic  harmonic  polynomial  solutions.  The 
regions  'R^L(QE]^',f,Th)  and  the  zero-contours  C°^(Qf  ;/;f,Tfc),  i  =  1,  2,  for 

k  =  1,. . .,5,  in  the  dement  f  for  the  five  3  x  3-cell  meshes  of  Fig.  19.  (a)  Mesh 
Tx ;  (b)  Mesh  Ta ;  (c)  Mesh  Ta ;  (d)  Mesh  T4 ;  (e)  Mesh  T, . 

Fig.  28.  Robustness  of  the  q%-superconvergence  regions:  Effect  of  the  surround¬ 
ing  mesh;  Elements  of  order  p,  harmonic  polynomial  solutions  of  degree  (p  -I- 1). 


sh;  Quadratic  dements,  cubic  harmonic  polynomial  solutions.  The 
(QEi  J;f,Tfc)  and  the  zero-contours  Cg_(Q?\I;f,Tk),  i  =  1,2,  for 

in  the  dement  f  for  the  five  3  x  3- cell  meshes  of  Fig.  19.  (a)  Mesh 


33 


The  common  region  lRj^(QH ;  /; f,  {T,- }£_,)  from  the  various  3  x  3- cell  meshes  for 
elements  of  degree  p;  (a)  p  =  2;  (b)  p  =  3. 

Fig.  24.  Robustness  of  the  t;%-super convergence  regions:  Effect  of  the  surround¬ 
ing  mesh;  The  five  5  x  5-cell  meshes  (with  the  element  f  shaded  gray)  formed  by 
surrounding  the  corresponding  3  x  3- cell  meshes  of  Fig.  19  by  one  layer  of  dements, 
(a)  Mesh  T';  (b)  Mesh  T';  (c)  Mesh  T';  (d)  Mesh  7”;  (e)  Mesh  T,'. 

Fig.  25.  Robustness  of  the  tj%- super  convergence  regions:  Effect  of  the  surround¬ 
ing  mesh;  Quadratic  dements,  cubic  harmonic  solutions.  The  regions 

jV'j^  and  the  zero-contours  (Of;  I;  f,  TL),  *  =  1,  2,  for  k  = 

1.. ...5,  in  the  dement  f  for  the  five  5  x  5-rdl  meshes  of  Fig.  24.  (a)  Mesh 
77;  (b)  Mesh  (c)  Mesh  T';  (d)  Mesh  T(;  (e)  Mesh  7?. 

Fig.  26.  Robustness  of  the  q%-superconvergence  regions:  Effect  of  the  surround¬ 
ing  mesh;  Cubic  dements,  quartic  harmonic  solutions.  The  regions 

ft  (0s i  f;  r,  TL)  *nd  the  zero-contours  (Q?;  I;  f ,  77),  *  =  1,  2,  for  k  = 

1. .  . . ,  5,  in  the  dement  f  for  the  five  5  x  5- cell  meshes  of  Fig.  24.  (a)  Mesh 
17;  (b)  Mesh  17;  (c)  Mesh  77;  (d)  Mesh  7?;  (e)  Mesh  T,'. 

Fig.  27.  Robustness  of  the  17%-super convergence  regions:  Effect  of  the  surround¬ 
ing  mesh;  Quartic  dements,  quintic  harmonic  solutions.  The  regions 
ft?|X  (QH-I-f,TL)  and  the  zero-contours  C°g  i  =  1,  2,  for  k  = 

1. . .,  ,5,  in  the  dement  f  for  the  five  5  x  5-cell  meshes  of  Fig.  24.  (a)  Mesh  T[\  (b) 
Mesh  77;  (c)  Mesh  T';  (d)  Mesh  T«';  (e)  Mesh  17. 

Fig.  28.  Robustness  of  the  q%-super convergence  regions:  Effect  of  the  surround¬ 
ing  mesh;  Elements  of  degree  p,  harmonic  polynomial  solutions  of  degree  (p  +  1). 
The  common  region  ft  J;  f,  {T/}i=1)  from  the  various  5  x  5-cell  meshes  for 

dements  of  degree  p:  (a)  p  =  2;  (b)  p  =  3. 

Fig.  20.  Robustness  of  the  17 %-super convergence  regions:  Effect  of  the  surround¬ 
ing  mesh;  Mesh  Tj'  formed  by  surrounding  the  5  x  5-cell  mesh  7^  of  Fig.  24e  with 
5  more  layers  of  uniform  dements.  The  dement  t  is  shown  shaded  gray. 

Fig.  SO.  Robustness  of  the  17%-superconvergence  regions:  Effect  of  the  surround¬ 
ing  mesh;  Elements  of  degree  p,  harmonic  polynomial  solutions  of  degree  (p  +  1). 
The  regions  (gg;  J;f,TgQ  and  the  zero-contours  C°Z  t  =  1,2, 

for  the  mesh  shown  in  Fig.  29.  (a)  p  =  2;  (b)  p  =  3;  (c)  p  =  4. 
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Fig.  SI.  Robustness  of  the  17%- super  convergence  regions:  Effect  of  the  surround¬ 
ing  mesh;  Mesh  formed  T*  by  surrounding  the  5  x  5- cell  mesh  T,  of  Fig.  24e  with 
5  more  layers  of  distorted  dements.  The  dement  r  is  shown  shaded  gray. 

Fig.  32.  Robustness  of  the  17%-superconvergence  regions:  Effect  of  the  surround¬ 
ing  mesh;  Elements  of  degree  p,  harmonic  polynomial  solutions  of  degree  (p  +  1). 
The  regions  j&Ji  and  the  zero- contours  (Qf  *  =  1,2, 

for  the  mesh  shown  in  Fig.  31.  (a)  p  =  2;  (b)  p  =  3;  (c)  p  =  4. 

Fig.  S3.  17%-superconvergence  for  general  solutions  and  meshes:  An  arbitrary 
mesh  Th  with  a  5  x  5-cdl  mesh  of  Union- Jack  pattern  embedded  in  it. 

(a)  Union-Jack  cell  (shaded  gray)  in  the  center  of  the  5x5  periodic  submesh; 

(b)  Union-Jack  cell  (shaded  gray)  off-center  of  the  5x5  periodic  submesh; 

(c)  Union-Jack  cell  (shaded  gray)  at  the  boundary  of  the  5x5  periodic  submesh. 

Fig.  34.  T/%-supcr convergence  for  general  solutions  and  meshes;  Cubic  dements, 
quartic  harmonic  solution.  The  regions  jVQ  (Qh\I\t,Tk)  and  the  zero-contours 

r,Th),  i  =  1,2.  (a)  Union-Jack  cell  shown  in  Fig.  33a;  (b)  Union-Jack 

cell  shown  in  Fig.  33b;  (c)  Union-Jack  cell  shown  in  Fig.  33c. 

Fig.  35.  17%-superconvergence  for  general  solutions  and  meshes:  The  grid  of 
Fig.  2  with  the  dements  r  (shaded  gray)  and  the  mesh-cell  w*.  (a)  Mesh-cell  1; 
(b)  Mesh-cell  2;  (c)  Mesh- cell  3. 

Fig.  36.  17 %-8uper convergence  for  general  solutions  and  meshes:  Periodic  cells 

endosing  the  mesh- cells  of  Fig.  35.  (a)  Periodic  cell  enclosing  mesh-cell  1;  (b)  Pe¬ 
riodic  cell  endosing  mesh-cell  2;  (c)  Periodic  cell  endosing  mesh-cell  3. 

Fig.  37.  17%-superconvergence  for  general  solutions  and  meshes:  Neumann  bound¬ 
ary  value  problem  over  the  domain  and  the  actual  mesh  of  Fig.  1;  linear  dements, 
harmonic  solutions  with  point  singularities  outside  the  domam  (see  Section  5.6). 
The  region  7i3g?(U;I; r,Th)  for  the  Zj-derivative  of  the  solution  for  (a)  Elements 

r  in  mesh-cell  1  (Fig.  35a);  (b)  Elements  r  in  mesh-cell  2  (Fig.  35b);  (c)  Elements 
r  in  mesh-cell  3  (Fig.  35c). 

Fig.  38.  17%-superconvergence  for  general  solutions  and  meshes:  Periodic  bound¬ 
ary  value  problem  over  the  periodic  cells  of  Fig.  36;  linear  dements,  quadratic 
harmonic  polynomial  solutions  obtained  from  a  local  Taylor-series  expansion  of 
the  solutions  with  point  singularities  outside  the  domain  (see  Section  5.6).  The 
region  7&jP(QB;I’,r,T)  for:  (a)  Elements  t  in  periodic-cdl  of  Fig.  36a;  (b)  El¬ 
ements  r  in  the  periodic-cell  of  Fig.  36b;  (c)  Elements  r  in  the  periodic-cdl  of 
Fig.  36c. 
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Fig.  SO.  q%-super  convergence  for  general  solutions  and  meshes:  Neumann  bound¬ 
ary  value  problem  over  the  actual  mesh  of  Fig.  1;  quadratic  dements,  harmonic 
solutions  with  point  singularities  outside  the  domain  (see  Section  5.6).  The  region 
for:  (a)  Elements  r  in  mesh  cell  1  (Fig.  35a);  (b)  Elements  r  in 

mesh  cell  2  (Fig.  35b);  (c)  Elements  r  in  mesh  cell  3  (Fig.  35c). 

Fig.  40.  q%-superconvergence  for  general  solutions  and  meshes:  Periodic  bound¬ 
ary  value  problem  over  the  periodic  cells  of  Fig.  36;  quadratic  dements,  cubic 
harmonic  polynomial  solutions  obtained  from  a  local  Taylor-series  expansion  of 
the  solutions  with  point  singularities  outside  the  domain  (see  Section  5.6).  The 
region  for:  (a)  Elements  r  in  periodic  cell  of  Fig.  36a;  (b)  El¬ 

ements  t  in  the  periodic  cell  of  Fig.  36b;  (c)  Elements  r  in  the  periodic  cell  of 
Fig.  36c. 

Fig.  41.  Robustness  of  the  q%-superconvergence  points  in  practical  computa¬ 
tions:  The  mesh  of  Figure  1  with  a  mesh-cell  u/jf  (shown  shaded  gray)  and  the 
corresponding  mesh-cell  (enclosed  by  the  thick-line  perigram). 

Fig.  42.  Robustness  of  the  17%-superconvergence  points  in  practical  computations: 
The  mesh-cell  u/£  with  the  regions  iV’Jj  (Cg;  J;tT),  for  the  levels  17  =  10,  20,  30, 

for  dements  of  degree  p.  (a)  p  =  1;  (b)  p  =  2;  (c)  p  =  3. 

Fig.  43.  Robustness  of  the  77 %-super convergence  points  in  practical  computations: 
Random  sets  of  10%-Buper convergence  points.  The  mesh- cell  u>%  with  the  regions 
wgiQ*-,  I\tT)  and  the  sampling  points  for  dements  of  degree  p.  (a)  p  =  1;  (b) 

p  =  2. 
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Regular,  Chevron,  Union-Jack  patterns 

Superconvergence  points  for  the  class  of  harmonic  solutions 

The  vertices  of  the  triangle 

are  at  (0.0, 0.0),  (1.0, 0.0),  (1.0, 1.0) 

Points 

P=  1 

p  =  2 

P- 3 

p  =  4 

p  =  5 

1 

0.5000000000, 

0.0000000000 

0.2113248554, 

0.0000000000 

0.0000000000, 

0.0000000000 

0.0118976218, 

0.0000000000 

0.0000000000, 

0.0000000000 

2 

0.7886751346, 

0.0000000000 

0.3752181980, 

0.0858825088 

0.2821230756, 

0.0000000000 

0.1726731646, 

0.0000000000 

3 

0.7500000000, 

0.0458758546 

0.5000000000, 

0.0000000000 

0.7178769244, 

0.0000000000 

0.5000000000, 

0.0000000000 

4 

0.7500000000, 

0.4541241452 

0.4715185566, 

0.3241733208 

0.9881023782, 

0.0000000000 

0.8273268354, 

0.0000000000 

5 

1.0000000000, 

0.0000000000 

0.1195044522, 

0.0224992688 

1.0000000000, 

0.0000000000 

6 

0.8638176260, 

0.6824066036 

0.7053452590, 

0.0071869378 

0.4640525316, 

0.0103019842 

7 

0.3095295948, 

0.2201118672 

0.2328766474, 

0.15749286240 

8 

0.6402277154, 

0.5196142358 

0.4734766188, 

0.3874100850 

9 

0.9063362308, 

0.7871241952 

0.7356302926, 

0.6425966766 

10 

0.9343659198, 

0.7080184052 

Table  la.  Classical  superconvergence  points  for  periodic  meshes:  Elements  of 
degree  p,  class  of  harmonic  solutions.  The  superconvergence  points  for  the  xa- 
derivative  of  the  error  and  the  class  of  harmonic  solutions  for  p  =  1,  2,  3,  4,  5 
for  the  element  r  (shaded  gray  in  Fig.  4)  for  the  Regular,  Chevron,  Union-Jack 
patterns. 
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Criss-Cross  pattern 

Superconvergence  points  for  the  class  of  harmonic  solutions 

The  vertices  of  the  triangle 

are  at  (0.0, 0.0),  (1.0, 0.0),  (0.5,  0.5) 

Points 

P-  1 

warn 

asm 

P  =  4 

mam 

1 

0.5000000000, 

0.0000000000 

0.2113248654, 

0.0000000000 

0.1464466094, 

0.0000000000 

0.0761088279, 

0.0000000000 

0.0488705882, 

0.0000000000 

■ 

0.7886751346, 

0.0000000000 

0.5000000000, 

0.0000000000 

0.3397707358, 

0.0000000000 

0.2355765633, 

0.0000000000 

B 

0.8535533906, 

0.0000000000 

0.6602292641, 

0.0000000000 

0.5000000000, 

0.0000000000 

0.5000000000, 

0.5000000000 

0.9238911721, 

0.0000000000 

0.7844234366, 

0.0000000000 

0.5000000000, 

0.4887686394 

0.9511294118, 

0.0000000000 

6 

0.5000000000, 

0.0000000000 

Table  lb.  Classical  superconvergence  points  for  periodic  meshes:  Elements  of 
degree  p,  class  of  harmonic  solutions.  The  superconvergence  points  for  the  xx- 
derivative  of  the  error  and  the  class  of  harmonic  solutions  for  p  =  1,2,  3,  4,  5  for 
the  element  Tj  for  the  Criss-Cross  pattern. 
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Criss-Cross  pattern 

Superconvergence  points  for  the  class  of  harmonic  solutions 


The  vertices  of  the  triangle  are  at  (0.0, 0.0),  (1.0, 0.0),  (0.5,  0.5) 


Points 

P  =  1 

P  =  2 

msm 

P  =  4 

P  =  5 

1 

0.5000000000, 

0.5000000000 

0.7886751346, 

0.5000000000 

0.8750000000, 

0.2834936534 

0.7621420044, 

0.3163540210 

0.8368698019, 

0.2268165908 

2 

0.8750000000, 

0.5000000000 

0.9145863153, 

0.3966634454 

0.7540881642, 

0.3397926069 

3 

0.8750000000, 

0.7165063710 

0.8970792676, 

0.5000000000 

0.6883014452, 

0.3781734456 

4 

0.5000000000, 

0.5000000000 

0.9145863153, 

0.6033365546 

0.7811733816, 

0.4582804499 

5 

0.7621420044, 

0.6836459789 

0.9334689552, 

0.5000000000 

6 

0.7811733816, 

0.5417195501 

7 

0.6883014452, 

0.6218265544 

8 

0.7540881642, 

0.6602073930 

9 

0.8368698019, 

0.7731834092 

10 

0.5000000000, 

0.5000000000 

11 

0.7760422930, 

0.5000000000 

Table  lc.  Classical  superconvergence  points  for  periodic  meshes:  Elements  of 
degree  p,  class  of  harmonic  solutions.  The  superconvergence  points  for  the  xx- 
derivative  of  the  error  and  the  class  of  harmonic  solutions  for  p  =  l,  2,  3,  4,  5  for 
the  element  t2  for  the  Criss-Cross  pattern. 


Superconvergence  points  for  the 
class  of  general  solutions 

Regular  pattern 

Points  from 
Table  la 

a 

a 

B 

D 

P- 5 

1 

yes 

yes 

no 

no 

no 

2 

yes 

no 

no 

no 

3 

no 

yes 

no 

no 

4 

no 

no 

no 

no 

5 

no 

no 

no 

6 

no 

no 

no 

7 

no 

no 

8 

no 

no 

9 

no 

no 

10 

no 

Table  2a.  Classical  superconvergence  points  for  periodic  meshes:  Elements  of 
degree  p,  1  <  p  <  5,  class  of  general  solutions.  This  Table  indicates  which  of 
the  superconvergence  points  for  the  class  of  harmonic  solutions  from  Table  la  are 
also  superconvergence  points  for  the  class  of  general  solutions  in  element  r  in  the 
Regular  pattern. 
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Superconvergence  points  for  the 
class  of  general  solutions 

Chevron  pattern 

Point*  from 
Table  la 

P=  1 

P=  2 

CO 

II 

a. 

•e 

II 

p  =  5 

1 

ye* 

no 

no 

no 

no 

2 

no 

no 

no 

no 

3 

no 

yes 

no 

no 

4 

no 

no 

no 

no 

5 

no 

no 

no 

6 

no 

no 

no 

7 

no 

no 

8 

no 

no 

9 

no 

no 

10 

no 

Table  2b.  Classical  superconvergence  points  for  periodic  meshes:  Elements  of 
degree  p,  1  >  p  >  5,  class  of  general  solutions.  This  Table  indicates  which  of  the 
superconvergence  points  for  the  harmonic  solutions  from  Table  la  are  also  super¬ 
convergence  points  for  the  class  of  general  solutions  in  element  r  in  the  Chevron 
pattern. 
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Superconvergence  points  for  the 
class  of  general  solutions 

Union- Jack  pattern 

Points  from 
Table  la 

P=  1 

P  =  2 

P  =  3 

P  =  4 

P=  5 

1 

no 

no 

no 

no 

no 

2 

no 

no 

no 

no 

3 

no 

no 

no 

no 

4 

no 

no 

no 

no 

5 

no 

no 

no 

6 

no 

no 

no 

7 

no 

no 

8 

no 

no 

9 

no 

no 

10 

no 

Table  2c.  Classical  superconvergence  points  for  periodic  meshes:  Elements  of 
degree  p,  1  <  p  <  5,  class  of  general  solutions.  This  Table  indicates  which  of  the 
superconvergence  points  for  the  harmonic  solutions  from  Table  1  axe  also  supercon¬ 
vergence  points  for  the  class  of  general  solutions  in  element  tx  in  the  Union-Jack 
pattern. 
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Superconvergence  points  for  the 
class  of  general  solutions 

Cris8*Cross  pattern 

Points  from 

Table  lb 

(For  the  element  rx) 

P  =  1 

P  =  2 

P  =  3 

p  =  4 

p  =  5 

1 

yes 

yes 

no 

no 

no 

2 

yes 

yes 

no 

no 

3 

no 

no 

no 

4 

no 

no 

no 

5 

no 

no 

6 

no 

Table  2d.  Classical  superconvergence  points  for  periodic  meshes:  Elements  of 
degree  p,  1  <  p  <  5,  class  of  general  solutions.  This  Table  indicates  which  of 
the  superconvergence  points  for  the  class  harmonic  solutions  from  Table  lb  are 
also  superconvergence  points  for  the  class  of  general  solutions  in  element  tx  in  the 
Criss-Cross  pattern. 
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Superconvergence  points  for  the 
class  of  general  solutions 

Criss-Cross  pattern 

Points  from 

Table  lc 

(For  the  element  r3) 

P=  1 

H 

p  =  3 

B 

H 

1 

no 

no 

no 

no 

no 

2 

no 

no 

no 

3 

no 

no 

no 

4 

no 

no 

no 

5 

no 

no 

6 

no 

7 

no 

8 

no 

9 

no 

10 

no 

11 

no 

Table  2e.  Classical  superconvergence  points  for  periodic  meshes:  Elements  of 
degree  p,  1  <  p  <  5,  class  of  general  solutions.  This  Table  indicates  which  of 
the  superconvergence  points  for  the  harmonic  solutions  from  Table  lc  are  also 
superconvergence  points  for  the  class  of  general  solutions  in  element  r3  in  the 
Criss-Cross  pattern. 
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Robustness  with  respect  to  the  type  of  solution 


*%(*) 

,l,f) 

Pattern 

Regular 

Chevron 

Union-Jack 

Criss-Cross 
(element  Tj ) 

Criss-Cross 
(element  ra) 

Points  from 
Tables  la,  lb,  lc 

rj % 

f}% 

if% 

1 

00 

100 

27 

00 

18 

2 

100 

27 

00 

3 

21 

83 

45 

4 

25 

53 

100 

Table  3.  Robustness  of  the  superconvergence  points  with  respect  to  the  type  of 
solution:  Quadratic  elements,  cubic  general  polynomial  solutions.  Value  of  tj%- 
error  in  the  Xj -derivative  of  the  solution,  at  the  various  points  corresponding  to 
the  superconvergence  points  for  harmonic  solutions. 
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Robustness  with  respect  to  the  distortion  of  the  mesh 


:»  mjx  0(4;  /;  C,  1,  f), 

f  m  F(i),  i'  =  f(m) 


Mesh  1  (Fig.  14a) 

Mesh  2  (Fig.  14c) 

Mesh  3  (Fig.  14d) 

Point*  from 
Ifcblela 

Harmonic 

Solution 

General 

Solution 

Harmonic 

Solution 

General 

Solution 

Harmonic 

Solution 

General 

Solution 

«J% 

9% 

iflf, 

6% 

3% 

1 

18.1 

28.0 

8.7 

75.2 

10.7 

94.1 

2 

19.2 

45.7 

16.8 

61.7 

37.1 

55.5 

3 

14.1 

66.8 

12.2 

50.7 

26.8 

48.8 

4 

2.7 

17.1 

9.6 

61.5 

22.9 

36.9 

Table  4.  Robustness  of  the  superconvergence  points  with  respect  to  the  distortion 
of  the  mesh:  Quadratic  elements,  cubic  polynomial  solutions.  For  the  meshes  of 
Fig.  13a,  Fig.  13c,  Fig.  13d,  obtained  by  displacing  the  central  node  of  the  Regular 
pattern,  the  value  of  17%-error  for  the  xx -derivative  of  the  error  is  given  at  the 
points  in  the  element  r  which  correspond  to  the  superconvergence  points  for  the 
undistorted  Regular  pattern  and  the  class  of  harmonic  solutions,  for:  (a)  Harmonic 
solutions;  (b)  General  solutions. 
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9%-error  at  a  random  set  of  10%-superconvergence  points 

s  ©(*;  «,  h,  r)  •  100% 

Point 

Value  of  q%  for  u  =  Uj 

Value  of  9%  for  u  =  Uj 

4.422 

3.592 

2.027 

6.461 

8.512 

4.917 

2.463 

0.671 

5 

2.801 

3.449 

mm 

5.884 

12.513 

mm 

3.369 

6.129 

8 

2.862 

4.335 

9 

2.594 

1.559 

10 

2.527 

1.489 

11 

1.859 

1.914 

12 

2.337 

1.433 

13 

6.109 

5.971 

14 

10.337 

11.140 

15 

6.456 

5.545 

16 

4.648 

11.874 

17 

7.086 

2.353 

18 

9.935 

14.901 

19 

9.848 

7.306 

2.779 

0.875 

21 

2.882 

4.376 

22 

6.195 

5.591 

Table  5a.  Robustness  of  the  ij%-superconvergence  points  in  practical  compu¬ 
tations:  Relative  error  at  the  10%-superconvergence  points  for  the  Xj-derivative 
for  linear  elements  (p  =  1).  These  points  were  randomly  selected  from  the  10%- 
superconvergence  regions,  as  shown  in  Fig.  43a. 


J7%-error  at  a  random  set  of  10%-superconvergence  points 

» 1%  =  ©(*;  £7; K  T ) '  100% 

Point 

Value  of  i}%  for  u  =  Uj 

Value  of  17%  for  u  =  u3 

11.936 

2.963 

8.743 

0.094 

4.372 

6.474 

0.582 

5.564 

5 

4.585 

1.350 

6 

'  ,/f 

6.036 

7 

2.086 

6.719 

8 

1.944 

6.918 

9 

7.722 

0.882 

10 

4 1-4 |j,  .;y1 

6.854 

11 

5.764 

0.962 

12 

5.331 

0.705 

13 

4.308 

0.852 

14 

1.447 

7.988 

15 

6.811 

2.824 

16 

6.317 

3.555 

17 

6.956 

0.047 

18 

8.101 

2.136 

19 

7.564 

4.765 

20 

6.637 

0.570 

21 

3.569 

1.081 

22 

5.423 

1.578 

23 

8.373 

0.432 

24 

7.065 

5.188 

25 

1.589 

5.472 

26 

4.162 

1.309 

27 

7.696 

0.494 

Table  5b.  Robustness  of  the  r/%-superconvergence  points  in  practical  compu¬ 
tations:  Relative  error  at  the  10%-superconvergence  points  for  the  Xj-derivative 
for  quadratic  elements  (p  =  2).  These  points  were  randomly  selected  from  the 
10%-superconvergence  regions,  as  shown  in  Fig.  43b. 
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Fig.  2a 


Fig.2c 


Fig.  3b 


Fig.  4a 


Fig.  4c 


Harmonic  solution  P  =  1 
SIMPLIFIED  L2  APPROACH 


0.0% 

10.0% 

20.0% 

30.0% 

Function  Q1 
Function  Q2 


Fig.  5a 


^ .0.0% 
mio.o% 

gg-20.0% 

&50.0% 

Harmonic  solution  P=  2 
SIMPLIFIED  L2  APPROACH 


Function  Q1 
Function  02 


0.0% 

10.0% 

20.0% 

30.0% 


Function  Q1 
Function  Q2 

Fig.  5c 


Harmonic  solution  P  =  3 
SIMPLIFIED  L2  APPROACH 


to.0% 

10.0% 

20.0% 

30.0% 

Harmonic  solution  P=  4  _  Function  Q1 

SIMPLIFIED  L2  APPROACH  _  Function  Q2 

Fig.  5d 


Harmonic  solution  P=  5 
SIMPLIFIED  L2  APPROACH 


Fig.5e 


0.0% 

■ 10.0 % 

■20.0% 

■30.0% 

Function  Q1 
Function  Q2 


I 


0.0% 

10.0% 

20.0% 

30.0% 

Harmonic  solution  P=  1  _  Function  Q1 

SIMPLIFIED  L2  APPROACH  _  Function  Q2 

Fig.  6a 


Bom 

iom 

20M 

30M 


Function  Q1 
Function  Q2 

Fig.  6b 


Harmonic  solution  P  =  2 
SIMPLIFIED  L2  APPROACH 


Harmonic  solution  P  =  3 

SIMPLIFIED  L2  APPROACH 

Fig.  6c 


10.0% 

10.0% 

20.0% 

30.0% 

_  Function  Q1 

_  Function  Q2 


10.0% 

20.0% 

30.0% 


Harmonic  solution  P  =  4 

SIMPLIFIED  L2  APPROACH 


Fig.  6d 


Function  Q1 
Function  Q2 


Harmonic  solution  P  =  5 
SIMPLIFIED  L2  APPROACH 


Fig.6e 


10.0% 

10.0% 

20.0% 

30.0% 

_  Function  Q1 

_  Function  Q2 


Ro.o% 

10.0% 

20.0% 

30.0% 

Harmonic  solution  P=  3  _  Function  Q1 

SIMPLIFIED  L2  APPROACH  _  Function  Q2 


Fig.  7c 


—-0.0% 

mio.o% 

1*20.0% 

■-50.0% 


Harmonic  solution  P  =  4  _  Function  Q1 

SIMPLIFIED  L2  APPROACH  _  Function  Q2 


Fig.  7d 


Harmonic  solution  P  =5 

SIMPLIFIED  L2  APPROACH 


0.0% 

10.0% 

20.0% 

30.0% 

Function  Q1 
Function  Q2 


Fig.7e 


rn 


Harmonic  solution  P  =  2 

SIMPLIFIED  L2  APPROACH 


Fig.  8b 


0.0% 

10.0% 

20.0% 

30.0% 

Function  Q1 
Function  Q2 


Harmonic  solution  P  =  4 

SIMPLIFIED  L2  APPROACH 


Fig.  8d 


- 0.0 % 

’10.0% 

r20.0% 

-30.0% 

Function  Q1 
Function  Q2 


Harmonic  solution  P  =  5 

SIMPLIFIED  L2  APPROACH 


Fig.8e 


—-0.0% 

mio.o% 

■ -20.0% 
SS-50.0% 

_  Function  Q1 

_  Function  Q2 


^0.0% 

mio.o% 

§§•20.0% 

General  solution  P=  1  ml~30.0% 

DIRECT  OPTIMISATION 

Fig.  9a 


{0.0% 
10.0% 
20.0% 
30.0% 

DIRECT  OPTIMISATION 

Fig.  9b 


General  solution  P  =  3 

DIRECT  OPTIMISATION 

Fig.  9c 


0.0% 

10.0% 

20.0% 

30.0% 


^ -0.0% 

■•20.0% 

■20.0% 

Harmonic  solution  P=  2  fH-30.0% 

SIMPLIFIED  L2 APPROACH 

Fig.  14a 


^ -0.0% 

mio.o% 

m-20.0% 

Harmonic  solution  P=  2  mL30.0% 

SIMPLIFIED  L2  APPROACH 

Fig.  14b 


Harmonic  solution  P  =  2 
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The  Laboratory  for  Nvawrkal  Analysis  is  an  integral  part  of  the  Institute  for  Physical 
Science  and  Technology  of  the  University  of  Maryland,  under  the  general  administration  of  the 
Director,  Institute  for  Physical  Science  and  Technology.  It  has  the  following  goals: 

•  To  conduct  research  in  the  mathematical  theory  and  computational  implementation  of 
numerical  analysis  and  related  topics,  with  emphasis  on  the  numerical  treatment  of 
linear  and  nonlinear  differential  equations  and  problems  in  linear  and  nonlinear  algebra. 

•  To  help  bridge  gaps  between  computational  directions  in  engineering,  physics,  etc.,  and 
those  in  the  mathematical  community. 

•  To  provide  a  limited  consulting  service  in  all  areas  of  numerical  mathematics  to  the 
University  as  a  whole,  and  also  to  government  agencies  and  industries  in  the  State  of 
Maryland  and  the  Washington  Metropolitan  area. 

•  To  assist  with  the  education  of  numerical  analysts,  especially  at  the  postdoctoral  level, 
in  conjunction  with  the  Interdisciplinary  Applied  Mathematics  Program  and  the 
programs  of  the  Mathematics  and  Computer  Science  Departments.  This  includes  active 
collaboration  with  government  agencies  such  as  the  National  Institute  of  Standards  and 
Technology. 

•  To  be  an  international  center  of  study  and  research  for  foreign  students  in  numerical 
mathematics  who  are  supported  by  foreign  governments  or  exchange  agencies 
(Fulbright,  etc.). 

Further  information  may  be  obtained  from  Professor  I.  BabuSka, Chairman,  Laboratory  for 
Numerical  Analysis,  Institute  for  Physical  Science  and  Technology,  University  of  Maryland,  College 
Park,  Maryland  20742-2431. 


